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Abstract 

Neutrino decay in vacuum has often been considered as an alternative to neutrino 
oscillations. Because non-zero neutrino masses imply the possibility of both neutrino 
decay and neutrino oscillations, we present a model-independent formal treatment of 
these combined scenarios. For that, we show for the example of Majoron decay that 
in many cases decay products are observable and may even oscillate. Furthermore, 
we construct a minimal scenario in which we study the physical implications of 
neutrino oscillations with intermediate decays. 
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1 Introduction 

The most favorable alternative for fast neutrino decay is to introduce an effec- 
tive decay Lagrangian which couples the neutrino fields to a massless boson 
carrying lepton number, i.e., a complex scalar field or a Majoron field [1-4]. 
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One possibility for such a Lagrangian for the case of Majoron decay is0 



Ant = Y.Y.9ij^lL^i,LJ, 



(1) 



where the i/'s are Majorana mass eigenfields and J is a Majoron field. 

Neutrino decay has been studied as an alternative to neutrino oscillations, 
especially for atmospheric [5-8] and solar [5,9,10] neutrinos. So far, either de- 
cay only or some sequential combination of decay and oscillations has been 
considered {e.g., MS W- mediated solar neutrino decay [10]). The values of in- 
dividual Qij^s or combinations were restricted, for instance, from double beta 
decay, pion and kaon decays, or supernovae [11-14]. For pure neutrino decay 
without oscillations the constraints on the coupling constants have been used 
to derive conditions for other parameters. For example, for atmospheric neu- 
trinos Am^ > 0.73 eV^ was found in Ref. [6]. Since this is incompatible with 
the known (or usually assumed) mass hierarchy of active neutrinos, decay into 
a sterile neutrino was suggested. However, massive neutrinos allow neutrino 
decay as well as neutrino oscillations. In this paper, we present a combined 
model-independent treatment of neutrino decay and neutrino oscillations. We 
refer to two different cases: 

Invisible decay: Decay into neutrinos or other particles which are not ob- 
servable {e.g., sterile decoupled neutrino^). 

Visible decay: Decay into neutrinos which are, in principle, observable {e.g., 
active neutrinos). 

Our paper is organized as follows: In Sec. 2, we will introduce a Majoron de- 
cay model for illustration, in order to show two important properties: First, 
as already mentioned before {e.g., in Ref. [5]), we may be able to detect active 
decay products. Second, we will show that in such a model, decay products 
can, in principle, oscillate. In Sec. 3, we will develop a model-independent op- 
erator framework for the calculation of transition probabilities and describe 

^ We are not interested in Vi transitions described by gu ^ 0, since they 

are forbidden for kinematical reasons. Furthermore, because the gij^s are usually 
assumed to be small, only the first order processes in the Lagrangian C\^t will be 
relevant to our discussion. 

^ In this paper, sterile decoupled neutrinos refers to neutrinos which can be de- 
scribed by a block diagonal mixing matrix such as 



where ns is the number of sterile neutrinos, t/sxs the active neutrino mixing matrix, 
and Ungxns the sterile neutrino mixing matrix. 
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its properties. Furthermore, we will show that the phase relationships given 
by the S'-matrix approach of quantum field theory are satisfied. The transition 
probability for invisible decay as a generalization of the neutrino oscillation 
formula will be derived in Sec. 4. As an example, we will apply this to at- 
mospheric neutrinos. In the following section, Sec. 5, this formalism will be 
used for visible decay in the approximation of almost stable decay products. 
Again, we will make an application to atmospheric neutrinos. In Sec. 6, we will, 
for illustration, construct a minimal model in which decay products oscillate. 
Finally, in Sec. 7, we will present a summary as well as our conclusions. 



2 Majoron models 

Neutrino decay through emission of a massless Nambu-Goldstone boson is 
currently among the most favorable alternatives for neutrino decay scenarios. 
Therefore, we will use it as an introductory example for our framework. 

2.1 Dynamics of Majoron decay 

For quantitative estimates we will investigate decay by Majoron emission de- 
scribed by an interaction Lagrangian with Yukawa (scalar or pseudoscalar) 
couplings (c/., Ref. [15]), 



Here ui and U2 are assumed to be Majorana neutrinos. Since the weak inter- 
action couples only to left-handed neutrinos and right-handed antineutrinos, 
and Majorana particles are identical to their anti-particles, Majorana neutri- 
nos and Majorana antineutrinos can only be distinguished by their spin. The 
matrix elements in the observer's rest frame are given by [15] 



int 



(3) 
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The differential decay rate for the decay U2 — i^i is 

■\M{E,){\ (6) 



dr2_^l 17111712 1 , . ./T^ x|2 



dEi AtT E2\P2\ 

with the kinematics in spherical coordinates (2;-axis in the direction of prop- 
agation) 

E2- Ei^ |P2 - pi| = -^/pi + P2 - 2|pi||p2| cos 6*. (7) 
The condition cos 9 <1 imphes for the allowed energy range of Ei that 



For relativistic neutrinos, ie., m2 <C £^2, this reduces to 

E2 



<Ei< E2. (9) 



We immediately see that for x ~ 1, i.e., Am^ <^ m^, the energy shift by decay 
is very small. For an energy resolution 

A£;,es > ^2 (1 - ^) (10) 

the decay products are detected in the same energy bin as the original particles. 

Now, we can integrate the differential decay rate given in Eq. (6) over the 
range determined by Eq. (8). Using Eqs. (4), (5), and (6) as well as assuming 
relativistic neutrinos, i.e., m2 <S E2, we obtain for the total decay rates [15] 



^ 2, 2 1 

•52 U-2 + -lna;+-- — 
V2 X x^ 2x^ 



and 



2.2 Visibility of decay products 



Solving Eq. (7) for cos^ shows that cos^ as a function of Ei always has a 
minimum at -Ei^min = "mi+m^ ^^S" ■'■)■ ^^^^ constrain 9. For 
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[eV] 



Fig. 1. The function cos^ = cos9{Ei) plotted for the (nonrclativistic) sample data 
mi = 1 eV, 1712 = 2 eV, and E2 = 10 eV. Since cos 61 < 1, this function is only 
defined below the horizontal line. 




/ L-l 



Fig. 2. A neutrino beam pointing from a source S to the detection area D with 
decay at X after traveling a distance I along the baseline L. The decay products of 
the beam cover the area F at the surface of the detector by a beam-spread of the 
angle 26. 

Table 1 

Typical energies (£^2), distances (d), and areas (F) covered by decay beam-spread 
for different types of neutrinos. Here it is assumed that m = 0{eV) and x = 0{1). 



Neutrino type 


£;2/MeV 


d/m 


F/m2 


Solar 


0.1 ~ 10 


-10^1 


10^0 ~ lO^i 


Atmospheric 


10^ ~ 10^ 


10^ ~ 10^ 




Accelerator 


10^ ~ 10^ 


10^ ~ 10^ 


10-16 _ 10-12 


Reactor 


1 ~ 10 


10 ~ 10^ 


10-12 ^ 10-8 


Supernova 


~ 10 


-10^1 


~ IOI6 
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relativistic neutrinos, i.e., 1712 <^ £^2, we obtain 

Ora^ = ^ - 1) • (13) 

Let us assume mi and m2 to be comparable by orders of magnitude, i.e., 
X ~ 1, which is plausible for active neutrinos for the currently assumed mass 
squared differences. For d, the distance to the detector, the area F covered by a 
neutrino beam at the detector due to a beam-spread by decay is approximately 
given by 

for X — 0{1). Table 1 shows typical values for several cases, where neutrino 
beams and not point sources are assumed. For a beam-spread at the detector 
(area F) much smaller than a usual detection area (D), F <^ D, everything will 
be visible, i.e., in the cases of atmospheric, accelerator, and reactor neutrinos. 

Let us now introduce a function 77 which describes the fraction of the decay 
products that will arrive at the detector by geometry, i.e., does not escape 
from detection by changing direction. We will first of all consider the case of 
a neutrino beam, and discuss the case of a point source later. We introduce 
this function in a quite general form applicable to any decay model such as 
a model involving more than two decay products. Thus, rj is a function of 
intrinsic decay properties such as energy and mass of the original particle Ei, 
mi, and of the decay products Ej, mj, as well as geometric properties such as 
length of the baseline L, decay position /, and detection area perpendicular to 
the beam direction D (see Fig. 2), i.e., 

T] = r]{Ei,mi,Ej,mj,L,l,D) = r]ij{L,l,D). (15) 

Note that the intrinsic properties of the beam and the detector can be inde- 
pendently folded with the transition probability, since they are not directly 
affected by decay. 

Since we assume that we cannot detect particles other than the neutrino decay 
product, we have to integrate over all momenta and energies except for the 
angle 9 {of.. Fig. 2) and the energy of the decay product E'. In general, these 
two parameters are to be integrated over the relevant detection region and the 
energy bin: 



r]ij{L,l,D) 



ptot 
ij 





d Fjj 


\ 




dE'd cos 9 

1 


V?min (cos 0)15 


d COS edE' 



(16) 



Here [£^min, -^'max] IS the energy bin range and (cos 9) d is the cosine of the 
maximal angle, such that the decay product will still hit the detector. Since we 
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Fig. 3. Tiie fraction 7721 = decay products of a solar neutrino beam hitting 

the detector as a function of the travehng distance from the surface of the Sun I. 
The sample data used are E2 = 10 MeV, mi = 0.02 cV, ma = 0.07 eV (c/., Ref. [5]), 
L = 10^^ m, and D = 10^ m^. Note that r/fi'^ ^ 1 for the Sun, because of the flux 
profile. 

only have one independent parameter in Majoron decay, the double-differential 
decay rate ^^^^q^^i will reduce to ^ times a 5-distribution determined by 
Eq. (7). Assuming that the energy bin width Aii^bin and the angle 9d are very 
small, i.e., the differential decay rate is approximately constant within the 
target region and the energy bin range, we obtain for 77 



ij 



d cos 9dE' 



AEun A{cos9)d, (17) 

cos6»=l, E'=Ei,in 



where Eun is the mean energy of the target bin considered, AE'bin the energy 
bin width of that bin, and A{cos6)d the cosine range, i.e., A{cos9)r) = 1 — 
{cos9)d, of the angle covered by the detector. From geometry (Fig. 2) we can 
determine A{cos9)d for the decay position / and a spherical detection area D: 

27r(L - 0^ A(cos 9)d ^ D. (18) 



Table 1 shows that for solar neutrinos the spread of a beam at the detector is 
normally larger than the detection area. Figure 3 shows a numerical evalua- 
tion of 7]2i for solar neutrino sample data and a very large detection area for 
illustration. One can see that in all cases the 1/{L — /) ^-dependence dominates 
in the range of about 6 • 10^° m to 9 • 10^° m over the effects of the couphngs 
or helicities, so that we obtain the same qualitative behavior. However, the 
Sun is a radially symmetric (47r) source. It is not a point source, because the 
production area is large compared to the distance to the Earth. Equation (13) 
indicates that for decay the maximal angle of deviation of the direction of 
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flight (9 ~ 10~^) is much smaller than the angle under which the Sun's core, 
where the neutrinos are produced, is observed from the Earth {6 ^ 10^'^). 
Hence, equally distributed small deviations of the direction of flight will al- 
ready be washed out by the flux profile. Besides that, after decay, the profile 
will still remain radially symmetric. Thus, for the case of the Sun, the function 
7^21 describing a beam can only serve as an example to show the qualitative 
behavior, and we conclude that 77^"" ~ 1. 

Supernova neutrinos are different from solar neutrinos, since they come from 
a far-distant point source, and hence from a well-defined direction. One might 
argue that it is impossible to distinguish neutrinos produced by decay in the 
direction of flight or by decay in slightly different paths due to the radially 
symmetric ffux profile, since decay scatters the neutrinos back into the di- 
rection of the detector. However, these two cases may be distinguishable for 
certain cases or sets of parameters because of the time resolution {e.g., different 
traveling times or different oscillation lengths due to different path lengths). 
Hence, supernova neutrinos might be the only nontrivial example of partly 
visible Majoron decay, i.e., rj < 1. In general, any far-distant beam source will 
show similar effects for Majoron decay. 

Eventually, we conclude that Majoron decay products are in most cases visible 
to a good approximation, i.e., riij{L,l, D) ^ 1. In our calculations, we will 
generalize this to rjij{L,l,D) ^ r]ij{L,D). Hence, rj will be assumed not to 
depend on the distance I. Moreover, for supernova neutrinos it should be a 
good approximation to use a step function for r]ij{L, I, D) according to Fig. 3. 



2.3 Neutrino oscillations of decay products 

In this subsection, we demonstrate that neutrino decay products in Majoron 
decay may oscillate. Since neutrino oscillations are a consequence of phase 
coherence, several issues connected to this subject need to be addressed for a 
combined treatment of neutrino decay and neutrino oscillations. We will do 
that in detail in Subsec. 3.2. 

In a charged current weak interaction, it is normally assumed that the phase 
coherence within a superposition of states as an in state is destroyed. This 
means that such a process produces a superposition of states with random 
relative phase shifts, i.e., the states within the superposition cannot interfere. 
Hence, cross sections referring to the individual states have to be summed in- 
coherently, i.e., by summation of probabilities and not amplitudes {of.. Prop- 
erty 4 in Subsec. 3.2). This comes from kinematical features related to the 
mass differences of accompanying leptons [16], i.e., the coherence length is 
relatively small due to the large mass differences. 
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In a neutral current weak interaction, the out state can be a coherent superpo- 
sition of states. One example is decay, producing a coherent superposition 
of different neutrino and antineutrino mass eigenstates [16] 

Wz) = ^ + \l^2)W2) + \l^3)W3)) (19) 

by an interaction Lagrangian 




Majoron decay resembles a neutral current weak interaction in the sense of 
a Lagrangian similar to Eq. (20). It couples a superposition of different mass 
eigenstates to a single particle, so that all eigenstates and their wave packets, 
respectively, experience the same energy shift AE transferred to the massless 
boson. From the wave packet treatment of neutrino oscillations we know that 
within a certain coherence length, the mass differences of neutrinos do not 
destroy the coherence of the mass eigenstates, i.e., neutrino oscillations. Hence, 
since in Majoron decay the energies of the mass eigenstates are all shifted by 
the same amount AE in decay, the mass differences will not be able to destroy 
the coherence of out states in this case eithei[^. Thus, the decay products can 
oscillate. 

Using a Majoron-like decay interaction Lagrangian, we know that to first order 
of the ^-matrix expansion we can write 



dT{in — > out) oc 



{out\ J d^xY^^gijPjPiJ \in) 



(21) 



This expression implies that the interaction Lagrangian destroys the incom- 
ing superposition completely, and creates a superposition of outgoing mass 
eigenstates with coefficients determined by the coupling constants and the 
coefficients of the incoming superposition. The outgoing states have a fixed 
relative initial phase such as the outgoing mass eigenstates in Z° decay [16]. 
In order to observe neutrino oscillations of decay products, we need to have 
at least two different mass eigenstates in the out state superposition. We im- 
mediately see from Eq. (21) that this implies that at least two of the (^j/s, say 
gi-^j^ and gi^j^, have to be non-zero for the corresponding states \vi^ and \vi^ 
with non- vanishing in state amplitudes Vi^ and z/jj. 



° Later in this paper, a more detailed proof of this analogy applied to our framework 
will be given in Property 9 in Subsec. 3.2. 
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3 Operators and properties 

In this section, we will introduce effective decay operators for a model-indepen- 
dent combined treatment of neutrino decay and neutrino oscillations. Decay 
can be described by decay rates oc gfj, which are calculated from the cor- 
responding Feynman diagrams. We will use such an approach to effectively 
integrate neutrino decay into neutrino oscillation scenarios. For this we do not 
need details of the specific model as long as neutrino decay can be described 
by a generic effective Lagrangian destroying one neutrino mass eigenstate and 
creating another one, i.e., Ci^t oc gijVjVi. This also includes multi-particle 
decay into decay products other than neutrinos, where all momenta and ener- 
gies of the undetected decay products are to be integrated over in the decay 
ratesQ. We will show that our operator framework satisfies the usual proper- 
ties of the S'-matrix approach and implements the correct phase relations. In 
certain cases, we will demonstrate this for Majoron decay. 



3. 1 Operators 

Let us introduce decay and propagation operators in terms of creation and 
annihilation operators, d) and d, respectively, as well as calculation rules for 
the transition probabilities. The symbols introduced in the operators will be 
explained thereafter. In the next subsection, we will show how the operators 
satisfy the expected properties. We define three operators: 

Definition 1 (Disappearance operator) is the transition operator gen- 
erally known as the "decay operator" . Effectively, it returns the amplitude for 
an undecayed state remaining undecayed after traveling a distance I along the 
baseline L: 



Definition 2 (Appearance operator) is the differential transition op- 
erator, which destroys an in state and creates an out state in [/, / -|- dl] along 



For decay models involving more than one outgoing neutrino we could use a 
similar framework witfi slightly different properties. For Majoron-like models higher- 
order processes are assumed to be suppressed, because the coupling constants are 
considered to be small. 
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the baseline L, i.e., a new state "appears 



-^Al.L) = EE^/¥ ^mALJ,D)d]a,e^^. (23) 

i j ^ * 

The phase ^ is a random phase taking into account the phase shift by addi- 
tional (not measured) particles produced in the decay, such as Majorons (cf.. 
Property 7 in Subsec. 3.2). The probability density, which is the square of the 
amplitude, will have to be integrated overt. 

Definition 3 (Propagation operator) £ is the operator propagating a state 
a distance I along the baseline ZQ; 

£{l) = Y.e^Y>{-iEil)a\ai. (24) 



In these definitions, rjij^L, I, D) is tlie geometrical function introduced in Eq. (15) 
and defined in Eq. (16). It describes the fraction of decay products which will 
still arrive at the detector, i.e., is not redirected from the detector by decay. 
The decay rate is defined as aij = ^ for i — > j decay, and r°- refers to the rest 

frame lifetime for that decay channel. Time dilation by the factor 7^ = Ei/rrii 
implies that 

pobserver _ pO . -1 _ ^ _ '^jj (nr\ 

« j * I 

The rate aj = J2j o^ij is the overall decay rate of the state i, so that B = ^ 
is the branching ratio for i j decay. The factor of two in the exponent of 
"D- and the square root in V+ comes from the fact that, besides some phase 
factors, amplitudes effectively behave like square roots of particles. However, 
as we will see in the next subsection, the above operators together with the 
definitions of the transition probabilities also implement the correct phase 
relations. Let us now define how to calculate the transition probabilities. This 
is to be done by successive application of decay and propagation operators: 

Definition 4 (Calculation of transition probabilities) The transition prob- 
ability P{h' a ^/s)^ = Pa/B exactly n intermediate decays is for n = given 
by 

P'^, = \{iy,\S{L)V4L)W^)\\ (26) 

^ Here we are using the non-relativistic operator, which is different from the rel- 
ativistic one only by an overall phase factor. This phase factor will cancel in the 
calculation of transition probabilities. For justification of this approximation see 
Ref. [17]. 
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and for n > given by 



L L 

P^0= j ■■■ J KM^L - l)V^{L - I) 

ll=0 ln=ln — l 

n f N n 

X n T^+{k, L)£{h)V.{k) |z/„) |2 n dk. (27) 
i=i ^ ^ 1=1 

Here I = J^^^ik- The order of and £ is arbitrary, since these operators 
commute as we will see below, in Eq. (32). 

In many cases, the total transition probability can be calculated as 

ax 

Pc.^ = E P:^^ (28) 

n=0 

where A^max denotes the maximal number of decays possible, depending on the 
decay model. The individual terms in Eq. (28) add up only if the particles 
produced by a different number of decays cannot be distinguished, such as by 
energy resolution, energy threshold or signature (flavor or spin). Depending 
on the problem, the terms may have to be split appropriately in genera^. 

Finally, we can introduce two simplified definitions for invisible as well as 
stable or almost stable visible decay products: 

Definition 5 {Invisible decay) For invisible decay products we can, in 
principle, only observe the undecayed particles themselves. Hence, we do not 
admit appearance operators, i.e., 

Piinvisible jdO /'on^ 
a/3 - ^al3- {^^> 



Definition 6 (Approximations for visible decay) Considering maximal 
one decay or lifetimes longer than we can neglect repeated decay terms 
(n > 1) and approximate the total transition probability in Eq. (28) by 

P ^ pO I pi _ pinvisible , pappearance /on^i 
-T a/3 — -T a/3 "I" -T a/3 — -T a/3 "r ^af3 ) ) 

where p^pp^^^'^'^'^'^ describes the "appearance" of new particles by decay. From 
Eq. (27) we can read off the first correction term P^^ for visible decay in 

^ For instance, we may even have two different P^^ coming from Vi Uj — > and 
Vi — > Vj — > Vk described by tlie appropriate decay rates. Anottier example is energy 
binning, wliere one may liave to be split into fractions of different energy ranges 
to be calculated from the differential decay rate, such as Eq. (6) for Majoron decay. 
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addition to p^^visibie . 

L 

Pip- J \{MS{L-l)'D-{L-l)V^{lL)E{l)V_{l)\v^)\^dl. (31) 
3.2 Properties 

In this subsection, we want to make the most important properties of our 

operator framework more transparent. Since we deal with decay of states in- 
stead of decay of particles, we need to show that this describes particle decay 
correctly while preserving interference properties for states. 

Property 1 (Creation and annihilation operators) The effective creation 
and annihilation operators, and a, respectively, are consistent with the S- 
matrix approach. 

The operator Ni = a\ai in Eqs. (22) and (24) is as usual the occupation number 
operator, whereas the operator T^- = a^cij in Eq. (23) is an effective transition 
operator coming from inserting the fermionic neutrino field expansions into an 
equation like Eq. (21). Here i and j refer to all quantum numbers, such as spin 
and momentum. Since we assume that conceptual issues, such as detection in 
a different energy bin or decay into an antiparticle, are dealt with by splitting 
Eq. (28) appropriately (c/., Def. 4), we do not need to take care of different 
quantum numbers here. 

Property 2 (Commutation relations) The operators T>_ and£ satisfy the 
commutation relation 

[P_,£] = 0, (32) 

whereas in general 

[P_,P+]^0, [P+,£:]^0. (33) 

These relations can be shown by direct calculation of the commutators, using 
the fermionic anticommutation relations {a^, a]} = 5ij. 

One can also derive Eq. (32) by introducing the notion of a complex mass 
square, i.e., 

= ml — iai. 

For £(l) — X^jexp [—lEil^ a-^Oj, where E'j ~ p + we can read off from 
Eqs. (22) and (24) that 

©_(/)£(/) =£{l) =S{l)V_{l). 
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This means that the disappearance and propagation operators can be com- 
bined by the introduction of complex mass squares. Therefore, the commutator 
reduces to the commutator of two scalars, i.e., the exponentials of the real and 
imaginary parts of the mass square. 

Property 3 (Decay into different channels) Our description for the de- 
cay of states corresponds to a description for the decay of particles, which is 
given by the differential equation system 

Here Tij is the rest frame decay rate for the particle decay i j and Ni is the 
number of particles i. From this we obtain for the disappearance of a particle 
Ni^o = 1, Njfi = for all j ^ i, and Vi = 'Ejj^i'^ij 

N,{t) = e-^^\ (35) 

and for the appearance of a stable decay product Ni^ = 0, Nj^ = 1 for all 
j 7^ i, and Fj = 

N.{t)= E ^fi (l - e-^^*) . (36) 



Neglecting propagation operatorsp"], we obtain from Eq. (26) for the disap- 
pearance probability of a mass eigenstate 

Pu = \{l^^\V^{l)W^)f = e-^\ (37) 

which is identical to the disappearance probability in Eq. (35) if Fj = 
which we have in the observer's rest frame. 

Neglecting propagation operators, we obtain from Eq. (31) for the appearance 
of a stable decay product {aj = 0) for exactly one decay 

P^J= r \{i^j\'D+il,L)V_il)\u,)fdl = v^(l-e-^'^), (38) 

Jl=0 0!i \ J 

which is identical to the appearance probability in Eq. (36) if 77 = 1 (everything 
detected) and Fj = which we again have in the observer's rest frame. 

In general, Eqs. (27) and (34) obey the same structure, since transitions among 
states are described by differential rates as in a Markov chain model. 

In this case, we are not interested in neutrino oscillation probabilities and phase 
factors will cancel anyway when squaring the probability amplitudes. 
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Property 4 (Coherent and incoherent summation of amplitudes) In 

quantum field theory, coherent and incoherent summation are often distin- 
guished in the calculation of Feynman diagrams [18]. The notion coherent 
refers to the summation of amplitudes 



dV oc 



and incoherent to the summation of squares of amplitudes, i.e., summation 
of probabilities 

dro^J2\Ai\\ 

i 

Basically, coherent summation applies to processes which are, in principle, 
indistinguishable and incoherent summation to processes which are, in prin- 
ciple, distinguishable. Our formalism implements the proper application 0/ co- 
herent or incoherent summation. 

To show this, we spUt the statement into more specific parts in what follows. 

Property 5 (Incoherent summation and repeated decay) Repeated de- 
cay can be described by the incoherent summation in Eq. (28). 

The differential decay rate in terms of the ^'-matrix is 

dV{in — > out) oc \{out \S\ in)f , (39) 

where S can be expanded as usually as [19] 

^ .-n °? ^ 

= E ^ / ■ ■ ■ / ^'^1 ■ ■ ■ d'XnT {Ant(2^l) ■ ■ ■ Ant(^n)} ■ (40) 



Thus, the in and out states select the corresponding S'-matrix expansion terms, 
which will evaluate to non-zero results. For a different number of decays we 
obtain different decay products characterized by different quantum numbers. 
Therefore, different ^'-matrix expansion terms will be selected in Eq. (39). 
This means that Feynman diagrams of different orders cannot interfere, even 
if we do not distinguish them by the measured decay products (often only 
the neutrinos but no other particles, such as Majorons, are measured). Hence, 
decay implies the principal possibility to measure the decay products and, 
thus, in terms of quantum mechanics acts as a measurement. Since repeated 
decay corresponds to higher order Feynman diagrams in quantum field theory, 
incoherent summation properly describes it. 

Property 6 (On-shell propagation in between decays) It is shown in 
Ref. [20] that Feynman propagators for virtual particles reduce to on-shell 
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propagators for macroscopic distance^^. Thus, individual decays in repeated 
decay can be treated independently without interference. More precisely, Feyn- 
man diagrams of the same order do not produce interference terms, because 
they can (in principle) be distinguished due to the large separation of vertices. 

Equation (28) together with Eq. (27) imply that decays are treated separately 
with decay rates calculated as for independent processes. 

Property 7 (Phase coherence between in and out neutrinos) Since at 
least a third unmeasured particle, e.g., the Majoron, is involved, the overall 
phase coherence between in and out neutrinos is destroyed. 

The random phase in Eq. (23) takes care of this. Since it does not appear in 
and cancels in P^^, it is only relevant for repeated decays for n > 1. It also 
ensures that terms from different individual decays do not produce interference 
terms, i.e., Feynman diagrams of the same order are summed incoherently in 
this case (c/., Property 6). 

Property 8 (Decay model and causality) The decay model used in Eqs. (26) 
and (21) properly implements causality. 

The operator D_ gives the probability amplitude for a state remaining unde- 
cayed until the position /. The operator P+ gives the conditional probability 
amplitude for decay in the interval [/, Z + c?/] for an undecayed state. Successive 
alternating applications of P_ and 1)+ give the (conditional) probability am- 
plitude for the decision path. This probability needs to be integrated to obtain 
the average over all possible decay positions / in a causal order, which is taken 
into account by the integration limits, i.e., a particle cannot decay before 
it even exists. We are integrating probability amplitudes for different decay 
positions incoherently, taking into account that undecayed particles cannot 
interfere with decay products. 

Property 9 (Neutrino oscillations of decay products) The appearance 
operator T>j^ gives the proper relative phase in a superposition of mass eigen- 
states as an out state. Kinematics induced by different masses does not result 
in destruction of coherence among the superimposed states such as in charged 
current weak interactions. 

The decay rate is to first order of the S-matrix expansion as usually defined 
as [19] 

(41) 



In this reference, this was explicitely shown for the creation and detection pro- 
cesses. 



rfm 



out) = — 
' T 



{out\ / d^x Citi<t\in) 
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Inserting a Majoron-like interaction Lagrangian {e.g., Eq. (1)) with the cor- 
responding field expansions yields for the decay rate from one flavor lu^) — 
T,iU^i into another one = J2j U^j \uj) 



^af3 — 



T 



(42) 



and from one mass eigenstate into another one {uj 

1 . 2 



r = — 



rp I «J I 



(43) 



Here the -Fj/s refer to transition functions depending on the properties of the 
states i and j as well as the Lagrangian and the used field expansions. Without 
loss of generality, for Dirac field expansions of fermions and a Klein-Gordon 
field expansion of the boson, we obtain for Lagrangians of the Majoron decay 
type 



/ (fx Antl^'i) 



(44) 



where Nij are normalization factors depending on energy, mass, and volume 
and A and B are some operators depending on the Lagrangian, e.g., combi- 
nations of charge conjugation and chiral operators. By introducing 0jj as the 
phase factor of Fij = \Fij \ e*'^'^ determined by Eq. (44), we can read off the 
following relation from Eq. (43): 



pi't'ij 



p. . 



As we will show below, for relativistic neutrinos (f)ij 
of the indices i and j. 



I.e. 



(45) 



is independent 



Without loss of generality, in our framework the rest frame transition proba- 
bility between two fiavor eigenstates for rjij = 1 {i.e., everything is detected) 
is given by 



a/3 



\{up\V+\u^)\'T 



« J 



T. 



(46) 



This can be obtained from Eq. (27) for an infinitesimally short baseline {dl — 
T) by using P_ £ i^lto zeroth order in the expansion of the exponentials 
as well as Tij = Applying Eq. (45) and the expansions of the fiavor 
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eigenstates, we finally obtain from Eq. (46) 



ra/3 = 



a/3 



1 

T 



(47) 



This is equivalent to Eq. (42) if 0,^ = is independent of the indices i and 
j. In this case, the appearance operator is identical to the expression derived 
from the -S'-matrix expansion and the different masses do not destroy phase 
coherence in a superposition of states. 

It remains to be shown that (f)ij is independent of the indices i and j for 
relativistic neutrinos. Since 4>ij = argFjj, any real factor in Fij such as, for in- 
stance, a normalization factor, does not affect (pij. We conclude from Eq. (44) 
that for real coupling constants the i and j dependence of (pij can only be 
influenced by the exponentials and the spinors. By spatial integration, the 
exponentials are transformed into ^-distributions implying momentum con- 
servation and, hence, cannot affect the phase. Without loss of generality, the 
spinors can be written in the Dirac representation as [19] 



1 





I 




1 



A^2(p^ - iv^) 



(48) 



where N\ and 1^2 are (real) normalization factors. For relativistic neutrinos 
we can assume the three momentum p to be identical for all mass eigenstates 
and absorb the mass dependence of the four-momentum in a correction of the 
energy ~ Then, they do not affect in this approximation, because 

the phases of these spinors only depend on the common three-momentum. 
Hence, is independent of the indices i and j. 



4 Invisible decay 



In this section, we will treat invisible decay with our operator framework, i.e., 
decay into particles, such as sterile decoupled neutrinos, which are in principle 
unobservable. In fact, this is the most often assumed neutrino decay scenario. 
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Fig. 4. The survival probability as a function of a and L/E for cos^ 6*23 = 0.30 
and Am^2 = 3-3 • lO'^ eV^ 

4-1 Transition probability for invisible decay 



The transition probability for invisible decay is calculated according to Eqs. (26) 
and (29). Since the operator V_ only gives real factors, the calculation for rela- 
tivistic neutrinos is quite straightforward and similar to the one for the general 
neutrino oscillation formula. We use |z/q,) = J2i U^i Wi) and define the following 
quantities: 

j:f^U^,Ul^U;,Up,, (49) 
Am? L (ml - mf) L 

vtiiL rrijL ( rrii mj\ L , . L . , 

Tii = — - — + ^ — ~ — + ^ — = (ai + a^ — . (51) 

Here Ei c:^ E has been used in the decay terms, because for decay the zeroth 
order terms in E'j ~ » + do not cancel, and therefore, the higher order 
terms can be neg 

lectedQ For more details on this kind of approximations see 



For neutrino oscillations the zeroth order terms cancel due to multiplication with 
their complex conjugates in the calculation of the transition probability. 
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Ref. [17]. 

After some calculations, not to be shown here, we finally obtain 



invisible ^ ^ ^ J^^e-^-- -4 T.'^Jtf sin^ A.^e" 

" . ' i>j 



pure decay 



PCPC onscrvmg 



2EE^^^'^sin2A,,e-^-. (52) 



« 1 

i>j 



Pop violating 

In our definition, pure decay admits neutrino fiavor mixing, but no neutrino 
oscillations (Ajj = for all indices i and j). 



4-2 Example: Application to atmospheric neutrinos 

In Ref. [7], neutrino decay was considered as an alternative to atmospheric 
neutrino oscillationsp^. Assuming Am^i = 0, 6^12 = 6'i3 = 0, and the CP- 
phase 5 = 0, we obtain from Eq. (52) for decay of z/2 into a sterile decoupled 
neutrino with the rate a = 02 the same transition probabilities as in Ref. [7]: 
For pure decay (Amgg = 0) we have 

^/.M= (46-^+4)', (53) 
and for decay and oscillations 



^t^t^ = 43 + 46 + 20^3523 COS ( ) e . (54) 



~2E 



The authors of Ref. [7] fitted Super-Kamiokande data in the case of decay and 
oscillations with negligible Am^g finding a = and cos^ 6 = cos^ 9[ 



723 



In this paper, we use the standard parameterization for the active neutrino mixing 
matrix [21] 

^ C12C13 S12C13 sise"**^^ 

^^3x3 = -S12C23 - Ci2S23Sl3e*'^ C12C23 - Sl2S23Sl3e*'^ S23C13 

\ S12S23 - Cl2C23Sl3e*'^ -C12S23 - Sl2C23Sl3e*'^ C23C13 J 

where Sjj = sin 9ij and Cij = cos 9ij and where 6ij are the vacuum mixing angles. 
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0.30. Pure decay without neutrino flavor mixing, i.e., cos^ 6 = 1, cannot de- 
scribe the atmospheric neutrino data [8,22]. Figure 4 shows the survival prob- 
ability P^^ as a function of the decay rate a and the sensitivity L/E. One 
can easily see that the transition from pure oscillation to pure decay is, in 
principle, determined by the curve ^ = 0{1). 



5 Visible decay 

In this section, we will assume that the neutrino decay products are, in princi- 
ple, observable. We know that we may either detect them separately from the 
undecayed neutrinos by different signatures or energies, or indistinguishably 
to the undecayed particles with the same signature in the same energy birf^. 
We postulated in Def. 4 that these conceptual differences manifest themselves 
in a problem- dependent splitting of Eq. (28) into appropriate parts. 



5.1 Transition probability for visible decay 

We use Eqs. (30) and (31) to first approximation for p^pp'^'^'^^^'^'^ ^ ^ ^j^^ case 
of maximal one decay or sufficiently small decay rates. For the calculation of 
the transition probability we make the following assumptions: 

• Stable decay products or <C for all i. 

• The fraction riij{L,l,D) of decayed particles not escaping detection by a 
change of direction does not depend on the decay position I, i. e., riij{L, I, D) = 
riij{L, D) = rjij is a (real) constant in /. In fact, we observed that r]ij{L, I, D) ^ 
1 in most cases of Majoron decay (c/., Subsec. 2.2). 

• The masses and decay rates are non-degeneratep^. 

For Tij and Aij we use the definitions in Eqs. (50) and (51). Applying Eqs. (30) 
and (31) and assuming that we cannot distinguish decay products from unde- 
cayed particles, i.e., we add and P^^, we have for Pq,^ 



This strongly depends on the detector properties. For example, for solar neutrino 
decay into antineutrinos with detection in Borexino or Super-Kamiokande it was 
discussed in Ref. [5]. 

The calculation of integrals changes for degenerate values of the a^'s and mj's (see 
calculation in Appendix A). Nevertheless, we will calculate limits with degenerate 
values later on. 
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P ~ pO _i_ pi pinvisible i pappearance 

-^a^+ / \{MS{L-l)V_{L-l)V+{l,L)S{l)V_{l)\i.^)fdl. (55) 



The invisible decay term, p^^i^ibie given by Eq. (52), i.e., 



invisible ^ ^^5^ ja/3^-r., 

i 3 



(56) 



i>j 



i>j 



The appearance term, pappearance evaluates after some algebra (c/., Ap- 
pendix A) to 



pappearance _ V^V^ / 
i j k I 



X 



^(r,,-r,,)^ + 4(A,, + A, 



^i^m) [(r,7 - r,fc) (e-^- cos(2Afc,) - e-^^' cos(2A;,- 

- 2 (Ai,- + A^jfc) (e-^'^^ sin(2Afei) - e'^^' sin(2A,,))" 

- '^{Kff.i) [{T,i - r,fc) (e-^^'= sin(2AfcO - e'^^^ sin(2A,, 

+ 2 (Ai,- + A^fc) (e-^^'^ cos(2Afei) - e"^^' cos(2A,,)) 
where K^Jli = U*j^UpjUakU^i is a generalization of J^"^. 

5.^ Limiting cases 



(57) 



Let us now consider some special cases for the general expression for the 
appearance term in Eq. (57). 



No oscillations. If we ignore neutrino oscillations, i. e., for degenerate masses, 
destruction of coherence in decay, or extremely small Am^'s, we will obtain 
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Fig. 5. The different parts of the survival probabihty plotted as functions of 
the sensitivity L/E ior the parameter values c^g = 0.3, a = 0.001^^, 7723 = 1, and 
Amis = 3.3- lO-^^eV^ [7]. 



in the hmit Am^ — > 



i j k I 



aj + ai - Qj - Qfc 



)• 

(58) 



No interference. For random phases within superpositions of mass eigen- 
states, i.e., 



(59) 



with a, b, and c some coefficients and ^1 and ^2 random phases, the interference 
terms in Eq. (57) vanish. We can include this by introducing SikSji in the sum 
of Eq. (57) and we finally obtain 



pappearance 

i j 



OLj — OLi 



2i£ 

E 



(60) 



Only one possible decay channel. Assume that there is only one possible 
decay channel, i.e., = ai = a^Q and all other a^i — 0. Then, we obtain 
the simplified formula 



pappearance 



(61) 
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5.3 Example: Application to atmospheric neutrinos 



Here we will extend the example from Subsec. 4.2, where Eq. (54) gives the 
oscillation formula for decay of V2 into an invisible decay product such as a 
sterile neutrino. Let us now assume that z/2 decays into the active neutrino 
v^. According to Eq.(55), we arc now looking for pappcarancc addition to 
P^txjx^^^^^i which can be calculated from Eq. (61) in this limit. If we assume 
that both undecayed and decayed neutrinos are detected indistinguishably, 
we will obtain 



p jjinvisible _j_ ^appearance (62) 

with P^^'^^'^ in Eq. (54) and 

pappearance ^ ^^^^2_^^2_^ _ ^-^^ (53) 



Here a is the decay rate of decay and E is the energy of the undecayed 

neutrino. Figure 5 shows the different parts of the survival probability for the 
parameter values given in the figure caption. One can easily see the effect of 
the increase of the survival probability for ^ ^ -. 



6 An example for neutrino oscillations of decay products 



In this subsection, we will construct a minimal example for observing neutrino 
oscillations among decay products. 



6.1 The scenario 



Since we want to create a superposition of out states in decay, we have to 
assume that at least two of the 5'ij's, i.e., cti/s, in the interaction Lagrangian 
are non-zero. We define a mass hierarchy ma > m2 > m\ and only allow the 
decay chain 1/3 — > i/2 — > i^i, i.e., 0:32 = q;3 7^ 0, 0:21 = q;2 7^ 0, and all the 
other Qjjj's are equal to zero. In order to demonstrate the physical effects we 
choose suitable parameters. Hence, we define the Am^'s quite close to each 
other, i.e., Am^i = 5 ■ 10~^ eV^, Am|2 = 3.3 • 10~^ eV^, maximal mixing, 
i.e., = ki) + ^2) + l^'a), and the a's close to each other, i.e., 
a2 = 1^^^ and 0:3 = yiUoT^" Furthermore, we assume that r]ij ^ 1 for all 
i, j and the CP-phase 6 — 0. 
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6.2 Transition probabilities 



Let us now look at the disappearance probability of electron neutrinos Pee- 
Prom Eq. (52) we obtain the invisible decay disappearance probability as 



pinvisible 7ee i Tee^-«2-fe i T^ee^-aa-^ 



+ 2Ji"|e-"2^ I 1 



+ 2J-e-"^t (i-2sin2'^"'3i 



'33^ 



+ 2 J23C 



^2 +"3 

2 B 



V 4 



(64) 



The appearance term is obtained from Eq. (57). Since in this case Kf^^^i — 
Ue^UejUekUei — | — tcrms in the sum of Eq. (57) are equal for 

simultaneous i ^ k and j / index exchanges. In order to make the result 
physically more transparent, we recognize that only a limited number of the 
terms in the sum are non- vanishing and split the formula into three parts: 



pappearance papp,l , papp,2 , pint 

ee ee ee ee " 

Here Pee"^'^ refers to the production of new uis by decay, 

p:r' = ^2121 (i - e 

PJ'PP'^ to the production of new i/2's by decay, 

,,2 r^ee "3 



papp,. 



-'^^3232 



-"3 I 



„ L 



0(2 — 0(3 



(65) 



(66) 



(67) 



and Pee to an interference term describing neutrino oscillations among the 
decay products ui and 1/2 (here K^22i — -^2132 used). 



pint _ A Ty-ee 
-'ee "^-'^^3221 



2 ^2 



X <^ as 



al + (Amii - Ami2) 

-23.L. ( Amni L 
e 2 E cos ' 



2 E , 



a2+"3 L ( ^rnln 

-e- . -cos^- 



"2 L 

e 2 B sin 



2 ^, 
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Fig. 6. The different parts of the probabihty p^Y^^'^^^"^'^^ plotted as functions oiLjE 
for the scenario data. 
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Fig. 7. The different parts of the probabihty Pee plotted as functions oi L/E for 

the scenario data. 

6.3 Results and analysis 



Figure 6 shows the different parts of the probability pappcarance ^-\^^ scenario 
data. The appearance of z/2 as decay product is described by PJJfP'^. It grows 
in the beginning, because it is dominated by decay. Since 0:2 > «3, the 

V2 is decaying faster than produced. Therefore, PlY^"^ falls again for large ^. 
The appearance of is determined by Plf^'^. Since vx is defined to be stable, 
PggPP'^ grows until it reaches its asymptotic value lim-c_,.oo PlY^'^ — -^^2121 — \- 

Describing neutrino oscillations of decay products, ie., and v^^ the in- 
terference term basically oscillates as a beat with the two frequencies 
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Fig. 8 
rates 0:3 and a2- 

determined by Amgg ± Am|^. The Amgg-dependency comes from the phase 
propagation before decay and the Am|]^-dependency from the phase propaga- 
tion after decay. 

Figure 7 shows the different parts of the total survival probabihty Pee- The 
invisible decay probability (disappearance) is the ordinary oscillation proba- 
bility damped by decay. Note that p^ppearance pmvisibie ^^^^ -^^^ ^^^y sensibly 

added if the decay products cannot be distinguished from the original ones 
in the detection process. In this case, the detection rates will be enhanced 
by a fixed amount at large L/E. Neutrino oscillations of undecayed particles 
will also vanish in this region. For small L/E there is a correction due to U2 
abundance by decay and in addition due to the interference term. 

Neutrino oscillations of decay products, induced by the oscillations of the 
interference term, can become unobservable due to suppression by the mixing 
angles or the decay rates different by some orders of magnitude. The first point 
is quite obvious, since the interference term is in general proportional to 



^ijkl — ^ai^' 



UakUti, where k^l, ^ {k, I). (69) 



This requires that there are coefficients of at least three different states in- 
volved as a product, i.e., all of them need to be quite large. The second point 
is physically obvious, but harder to see in our equations. Since the interfer- 
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Fig. 9. The probability p^^^^^'^'^^^'^ for different limiting cases plotted as functions 
oi L/E for the model data. 

ence describing neutrino oscillations of decay products can only be observed 
if two states are produced simultaneously in decay, the couplings to both of 
them need to be quite large. Large couplings imply large decay rates for both 
transitions. If one of them is too large, the corresponding state will have de- 
cayed before oscillations can be observed. Hence, the decay rates need to be 
large enough and of the same order of magnitude in order to obtain reasonable 
results. A measure for the magnitude of the interference effects is the absolute 
maximum of |-Pee*l- ^^Z- 8) we show the result of a numerical evaluation of 
this function for different values of the decay rates as and a2- One can easily 
verify that for too small or too different decay rates, this function gives small 
values, i.e., we do not see any interference effects. 

6.4 Application to limiting cases 

In Subsec. 5.2, we studied the appearance probability in certain limits. Let 
us now apply these to our scenario to see what happens if certain terms are 
assumed to be negligible. 



No oscillations. In this hmit, all Am^'s are equal to zero (c/., Eq. (58)). 
Here we can simply use 




+ 



pint 

ee 



(70) 
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Fig. 10. Possible decay scenario for four neutrinos. 



No interference. This limit is described by Eq. (60), which is equivalent 
to neglecting the interference term. We obtain 

jjappeajrance papp,! _j_ papp,2 (71) 

Figure 9 shows the different limits of the appearance term compared to the 
original one. In all cases, we obtain similar values for small and large L/E, since 
for L/E (Am^)^^ the oscillation terms are ineffective and for L/E 
they are suppressed by decay. In the limit of no oscillations, there is no oscillat- 
ing behavior of the curve. Nevertheless, there are interference terms coming 
from the coherent coupling of the Lagrangian to different mass eigenstates. 
Hence, in this case we obtain an enhancement in the appearance probability. 
In the limit of no interference, we obtain a curve similar to the original one, 
but without oscillations, which are described by the interference term. 



7 Sumniciry and conclusions 



In this paper, we have introduced a combined treatment of neutrino decay and 
neutrino oscillations in a formal model-independent operator framework. We 
have shown that it can be easily applied to special cases such as atmospheric 
neutrinos. We have also seen that for Majoron decay active neutrinos as decay 
products have to be taken into account in the detection rates and may even 
be indistinguishable from the original undecayed particles by the detector. In 
addition, in Majoron-like decay models neutrino decay products may oscillate. 

We conclude that the visibility of decay products in Majoron decay models 
may crucially affect the detection rates depending on the properties of the 
detector. For other decay models only a certain fraction of the decay products 
may arrive at the detector due to kinematics. Moreover, neutrino oscillations 
of decay products are observable for relatively large neutrino flavor mixing 
only. 

One can show that for the combination of neutrino decay and neutrino os- 
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cillations, the dimension of the parameter space is equal to |n(n — 1), where 
^n{n — 1) parameters come from decay and n{n — 1) from oscillations. Here 
n denotes the number of neutrino flavors. Hence, the number of parameters is 
increased when neutrino decay is included. Thus, especially in multi-neutrino 
scenarios, effects not consistent with the conventional neutrino oscillation 
framework could be described by a combined neutrino decay and neutrino 
oscillation approach. 

As a realistic scenario, we could, for example, use a four-neutrino framework 
with a decoupled sterile neutrino i/s (cf-, Eq. (2)) and a mass hierarchy 777.3 > 
7712 > iTT-i > f^s (c/, Fig. 10). In such a scenario, we have up to six decay rates 
in addition to the usual mixing parameters. From other conditions such as 
the observability of supernova neutrinos, we then would need to constrain the 
lifetimes for individual i/j's {e.g., v\ almost stable). Since the maximal number 
of possible decays is A^max = 3 in this scenario, the exact transition probability 
Pq,/3 without lifetime constraints would be determined by P°^, P\^, Pa^i ^iid 
Pip in Eq. (28). 

Since non-zero neutrino masses imply both neutrino decay and neutrino oscil- 
lations, neutrino decay may become interesting not only as an alternative to 
neutrino oscillations but also as a part of the whole neutrino scenario. Even 
small (jfjj's may affect detection rates of long traveling neutrinos such as su- 
pernova or primordial ones. Flux measurements {e.g., of supernova neutrinos) 
cannot always be used to extract neutrino mass limits reliably, because of 
different path lengths of undecayed and decayed neutrinos. Thus, dispersion 
may be mimiced by decay. 

As a final conclusion, our operator formalism is a generalization for a combined 
treatment of neutrino decay and neutrino oscillations could be of interest to 
those, who are interested in analyzing neutrino data in a more general way. 
In addition, it may also be applied in other oscillation scenarios outside of 
neutrino physics, where coherence is not destroyed in decay. 
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A Calculation of the appearance term 



We are using Eq. (31) as an approximation for the appearance term: 



pappearance 



1=0 



dl 



A^pl dl, 



(A.l) 



1=0 



where Aap is the (differential) transition amplitude. Applying Eqs. (22), (23), 
(24), and |z/a) = \ we obtain for the transition amplitude]^ 




an , 



(A.2) 



In this equation, we already have evaluated the creation and annihilation op- 
erators for readability. Now we use the approximation for relativistic neutrinos 

2 2 

Ei c:^ p + ^ For neutrino oscillations the three- momentum p only 

gives a common phase factor, which will cancel. For neutrino decay, the expo- 

are real, and thus Ei 



nentials e and e are real, and thus E'j ~ p ~ ii^ in the lowest order 

nontrivial approximation. Using the relativistic approximations as well as the 



definition Am^^ 



then yields 



111 



(A.3) 



*7^i Z— independent 

For the appearance probability we obtain 



dependent 



L 1. 

pappearance^ / \A^,f dl = f A^.A^.rf/ = ^ EE E f^^.f^^.f^^^f^^Z 



Q/3 

ijkl 



Q,+Q, Am? - . A„2 



(A.4) 



Here we already neglect the random phase factor e*^ from the operator since it 
will cancel anyway by taking its absolute value squared in the transition probability. 
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where the rjij^s are assumed to be reaL For further evaluation we use the 
integral {a, b real) 



for a + bi 0. Comparison with Eq. (A.4) shows that the latter condition 
implies 

aj + ai — ai — ak + i (m^ — rrtf — ntf + ^ 0. (A. 6) 

Since i^j/Xk^lhy the summation rules, this can only happen for {i = 
k) A (j = /) A (ctj = aj) for non-degenerate m's. Here we assume the a's to be 
non-degenerate for simplicity, so that we can apply Eq. (A. 5). Using the same 
abbreviations as for invisible decay in Eqs. (50) and (51), we can continue 
evaluating Eq. (A.4) with Eq. (A.5) by identifying a = = Vji - Tik 

and b = = 2 (A,, + A,,) 



T a ■ -\-ai Am 

^appearance = ^ VW^V^ p^" ^"e"^^" 

i j k I ^ 

X - Tik - 2i (A,-, + Am) / („.+„,_„._„,+i(Am|,-Am2,))^ _ 

(r,,-r,fc)V4(A,, + A,,)' V 

i j fc / (i j7 - i ife) 4 {Aji + Aki) 

T-) T-iappearancc • i i ^appearance crv ^appearance tti j. i 

Because F^J^ is real, we have — ^jir^^ . box cp real, y 

and 2; complex, it can be shown that 



z cos ( 



'^[z'] sin< 



^{iye^'^) = -^{y) COS (j) - ^{y) sin < 
We can use this to find 



(A.8) 
(A.9) 



3? (X(c + di)e'^'^) = ^(K)c cos 2 A - $J (ir)d sin 2 A 

-S(X)ccos2A-K(X)dsin2A (A. 10) 

for c, real. In comparison with Eq. (A. 7), we identify c = Tji — Tik and 
d = 2(Ajj -|- Aik). Applying this formula and re-grouping in and '^K 
yields the final result in Eq. (57). 
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